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Heavy ion collisions at RHIC produce direct photons at low transverse momentum, pv from 
1-3 GeV/c, in excess of the p+p spectra scaled by the nuclear overlap factor, Taa- These low pv 
photons have a large azimuthal anisotropy, V 2 - Theoretical models, including hydrodynamic models, 
struggle to quantitatively reproduce the large low pt direct photon excess and V 2 in a self-consistent 
manner. This paper presents a description of the low pr photon flow as the result of increased 
photon production from soft-gluon mediated q-q interactions as the system becomes color-neutral. 
This production mechanism will generate photons that follow constituent quark number, Uq, scaling 
of V 2 with an Uq value of two for direct photons, comparisons of the published PHENIX direct 
photon and identified particle V 2 measurements finds that nq-scaling applied to the direct photon 
V 2 data prefers the value Uq = 1.8 and agrees with Uq = 2 within errors in most cases. The 0-20% 
and 20-40% Au-|-Au direct photon data are compared to a coalescence-like Monte Carlo simulation 
that calculates the direct photon V 2 while describing the shape of the direct photon pv spectra in a 
consistent manner. The simulation, while systematically low compared to the data, is in agreement 
with the Au-|-Au measurement at pv less than 3 GeV/c in both centrality bins. Furthermore, this 
production mechanism predicts that higher order flow harmonics, Vn, in direct photons will follow 
the modified Wq-scaling laws seen in identified hadron Vn with an Uq value of two. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 


Direct photons are all of the photons produced in a col¬ 
lision excluding the products of hadronic decays. They 
are emitted throughout the evolution of the heavy ion 
medium, and because they are color-neutral they do 
not experience subsequent interactions with the medium. 
As a result, their spectrum provides a time-integrated 
picture of photon emission. Direct photons have var¬ 
ious sources, including prompt photons generated by 
early hard parton interactions, photons produced in 
the pre-equilibrium stage, and thermal photons radi¬ 
ated from either the quark gluon plasma (QGP) or the 
hadron gas stage (HG). In Figure [U Feynman diagrams 
of prompt photon production mechanisms, quark-gluon 
Gompton scattering, quark-anti-quark annihilation, and 
bremsstrahlung radiation, are shown. Prompt photons 
are created in p+p collisions and dominate the yield at 
high Pt in heavy ion collisions. Prompt photon pro¬ 
duction rates can be calculated using perturbative QCD 
(pQCD); quark-gluon Gompton scattering and quark- 
anti-quark annihilation have production rates of order 
asoi and bremsstrahlung radiation has a rate of order 
Q:|a. QCD thermal photons have the same production 
diagrams, shown in Figure [2 but with the partons ther- 
malized in the medium. In thermal photon pQCD calcu¬ 
lations, bremsstrahlung radiation is of order asa and can 
exceed the production from the Compton scattering and 
annihilation processes. HG thermal photons have analo¬ 
gous production mechanisms to the Compton scattering 
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and annihilation processes only with pions and p-mesons 
interacting instead of quarks and gluons. However, the 
production rates for thermal photons and other direct 
photons sources are not well constrained particularly in 
the non-perturbative regime. This makes separating the 
contributions of direct photons at low and intermediate 
Pt difficult. 

The PHENIX experiment discovered a large di¬ 
rect photon excess at low pT, from 1-3 GeV/c, in 
y/Sjqjq = 200 GeV Au-I-Au collisions at RHIC relative to 
the yields of direct photons in collisions scaled by the 
nuclear overlap factor, Taa Islsl- Subsequent analyses 
found that these low pT photons, again from 1-3 GeV/c, 
have a large azimuthal anisotropy with respect to the 
collision’s event plane Q. Preliminary results from the 
ALICE experiment at the LHC SM^st similar behav¬ 
ior in 2.76 TeV Pb-|-Pb collisions [3,0. Hydrodynamic 
models are able to describe the direct photon yield with 
initial temperatures of 300-600 MeV and thermalization 
times between 0.15-0.5 fm/c 0. Reproducing the large 
measured azimuthal anisotropies, V 2 , at these early times 
has proven difficult for hydrodynamic models |7H9j . This 
is because the large azimuthal anisotropies generated by 
hydrodynamic pressure gradients need time to develop. 
To address this puzzle some theories introduce delayed 
QGP formation [^, new sources of photon production 
involving stron g rn agnetic fields [l^ and initial state 
Glasma effects [Q, while others consider increased con¬ 
tributions from the hadron gas stage due to baryon- 
baryon and meson-baryon interactions 0[ll. 

In this paper, the sources of identified hadron az¬ 
imuthal anisotropies are considered to understand the 
origin of the similarly-sized direct photon V 2 - At low pT, 
bulk expansion dominates the hadronic V 2 while at high 
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FIG. 1: Feynman diagrams of prompt photon production by a) 
lation, and c) bremsstrahlung radiation off of an outgoing quark 


quark-gluon Compton scattering, b) quark-anti-quark annilii- 


q 



FIG. 2: A Feynman diagram of the quark-anti-quark anni¬ 
hilation interaction with a medium gluon producing a direct 
photon. 


Pt, hadrons from jet fragmentation dominate. In the 
intermediate pr region, from 1-3 GeV/c, the measured 
baryon and meson V 2 values split, with baryons reaching 
higher values of V 2 at higher values of px [la| . When the 
baryon and meson V 2 values are scaled by their number 
of constituent quarks, rig, a uniform behavior between 
baryons and mesons is seen El- Coalescence models are 
able to reproduce quark number scaling by assuming that 
hadron production is dominated by the recombination 
of flowing partons. They assume that thermalized co¬ 
moving quarks of a given px will coalesce into mesons and 
baryons with rig-times the px and rig-times the V 2 where 
rig = 2 for mesons and Uq = 3 for baryons. In this frame¬ 
work, energy-momentum conservation is maintained by 
the mean-field interaction resulting in soft gluon interac¬ 
tions with the medium [l^ . 

Similar mean-field or soft gluon interactions could me¬ 
diate quark-anti-quark annihilation as the system moves 
toward color neutrality, resulting in a large increase in 
photon production. These interactions, a diagram is 
shown in Figure [21 would produce photons from partonic 
processes late in the system’s evolution when quarks are 
flowing. One consequence of this production is that these 


photons should reproduce constituent quark number scal¬ 
ing with the value rig = 2 for direct photons. Fur¬ 
thermore, this model provides a testable prediction that 
higher order flow harmonics, Vn, in direct photons should 
follow the rig-scaling laws seen in identified hadron Vn [3 
again with rig = 2 for direct photons. 

Section El determines the rig for direct photons that 
best reproduces the quark number scaling seen in the 
identified hadron V 2 by using a analysis of existing 
data i0- Section uni details a coalescence-like Monte 
Carlo calculation that combined with the T^A-scaled p+p 
component is compared to the measured direct photon px 
spectrum and V 2 distribution. A two-component model is 
assumed where the low px direct photon excess is primar¬ 
ily the result of quark-anti-quark annihilation mediated 
by mean-field or soft gluon interactions as the system 
becomes color neutral. 


II. THE rig-SCALING OF IDENTIFIED 
HADRON AND DIRECT PHOTON V 2 

The elliptic flow of identified hadrons displays con¬ 
stituent quark number scaling in the 1-3 GeV/c px re¬ 
gion [21I, [22. In the q-q annihilation picture of direct 
photon production, this rig-scaling behavior should ex¬ 
tend to the direct photons with rig = 2. This is be¬ 
cause the rig-scaled V 2 reflects the underlying anisotropy 
of the quarks and therefore is common for all hadrons 
and photons produced from these coalescing quarks. At 
high Px-, this rig-scaling may breakdown as contributions 
from hard processes begin to dominate in both the direct 
photon and identified hadron spectra. Figure |3| shows a 
comparison of the direct photon V 2 0 with the charged 
pion, kaon and proton V 2 [l^ in the 0-20% and 20-40% 
= 200 GeV Au-I-Au collisions. The rig-scaled V 2 
as a function of the rig-scaled px and KEx are also pre¬ 
sented assuming that the rig value for direct photons is 
two. The agreement between the scaled direct photon V 2 
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and the pion, kaon and proton data is impressive despite 
the large systematic error bars on the direct photon mea¬ 
surement. The scaled pions, kaons, protons and photons 
agree at low KEx/nq in both centralities. At KET/uq 
above 1.7 GeV, the direct photon’s scaled drops below 
the pion values. This deviation can be understood as the 
result of the increased photon production by initial hard 
processes Q . Of particular note is how the direct photon 
and proton V 2 /Uq track together as they deviate from the 
pion values in the 20-40% centrality bin. This suggests a 
similar transition to the high hard scattering region 
for the scaled protons and photons. While the 0-20% 
proton V 2 does not extend high enough in KEr/nq, pro¬ 
tons in the 0 - 20 % centrality are also expected to break 
Uq-scaling at high KEr/uq and deviations are seen in 
the 10 - 20 % bin 

A analysis is undertaken to determine if = 2 
best produces the agreement between the direct photon 
and the rig-scaled identified hadron V 2 data. This is done 
in two ways. In Section fll A[ the datasets are compared 
directly. In Section HlBi the rig-scaled identified hadron 
V 2 are fit and the direct photon V 2 are compared to that 
function. 


A. comparison between the direct photon and 
rig-scaled hadron data 


A x^ comparison is performed between the V 2 for direct 
photons to the rig-scaled hadron data. The compar¬ 
ison of the direct photon and identified hadron data is 
calculated according to 


X 
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E E E 

Cent. tt.K.p KEt jn^ 
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( 1 ) 


where U 27 is the direct photon V 2 -, V 2 h is the identified 
hadron V 2 for each of the summed hadrons, rr, K and p. 
The x^ is summed over the 0-20% and 20-40% centrali¬ 
ties comparing the rig-scaled pion, kaon and proton V 2 juq 
values to the direct photon V 2 /nq^ where rig-,, is the only 
parameter. Determining the photon and hadron uncer¬ 
tainties, (T-y and ah^ is complicated because the published 
systematic errors for both the direct photons and identi¬ 
fied hadrons combine both point-to-point and correlated 
systematic errors BUS- To address this the x^ analysis 
is performed in two ways. In one case, the quadrature 
sum of the statistical and systematic errors for direct 
photons and the identified hadron uncertainties is used, 
O’ = cTstat © cTsys- This assumes that the systematic er¬ 
rors are uncorrelated. Another x^ analysis assumes that 
the systematic errors are fully correlated and the photon 
and hadron uncertainties are limited to their statistical 
errors, a = (Tstat- In both cases, the comparison of a 
given pair of direct photon and hadron data points are 
included in the x^ calculation only if the KE^/nq values 
are within 0.1 GeV/c of each other. An example of this 
data comparison over the full range in REx/nq is shown 




pyriq (GeV/c) 


p^tiq (GeV/c) 



0.15 

0.10 

0.05 

0.00 


KEt/m (GeV) 


KEjIn (GeV) 


FIG. 3: (color online) The tt (blue triangles), K (open ma¬ 
genta squares), p (red circles) and direct photon (open green 
crosses) V 2 as a function of pr in central 0-20% (a) and mid¬ 
central 20-40% (b) Au-l-Au collisions at ~ 200 GeV. 

Panels (c) and (d) show the as functions of pxjuq for 

0-20% and 20-40% respectively. Panels (e) and (f) show the 
U2 scaled by the number of constituent quarks, rig, as a func¬ 
tion of KEt/ riq, again for 0-20% and 20-40% centralities. For 
direct photons, rig = 2 is assumed. In panels (a), (c) and (e), 
the 0-20% V 2 values are scaled by 1.6 for better comparison 
to the 20-40% results. Error bars and shaded boxes around 
points represent their statistical and systematic uncertainties 
respectively 0 , 0 - 


in Figure m where the photon-to-identified hadron data 
comparison plots with Uq^ = 2 are presented. A x^ of 
16.28 is calculated using the quadrature sum of the sta¬ 
tistical and systematic errors for the photon and hadron 
uncertainties with 35 degrees of freedom, NDF, and a 
reduced x^, /NDF, of 0.47 is found. As a result of re¬ 
quiring photon-hadron matching in KEp /iig, the number 
of degrees of freedom of the x^ calculation changes as rig-,, 
varies. This leads to a discontinuous x^ distribution as a 
function of rig.y, as seen in Figure [SJ 

Figure!^ (a) shows the x^ versus Uq^ when statistical 
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FIG. 4: (color online) Example plots of the input data used for the calculation of comparing the v^juq vs KET/uq for 
identified hadrons (red circles) [2^ and direct photons (black squares) [j using the quadratic sum of the statistical and sys¬ 
tematic errors. Here, nq-y = 2 is assumed for direct photons. The 0-20% (top row) and 20-40% (bottom row) y^sTTj^T = 200 GeV 
Au+Au results are shown. Pions (left column), kaons (middle column) and protons (right column) are separately plotted with 
the direct photon data over the full KEr/riq range. The data are included in the calculation only if the identified hadron 
and direct photon KET/riq values are within 0.1 GeV/c. The x^ is calculated using the variation between direct photon and 
identified hadron V 2 juq in all six plots. Error bars represent the statistical and systematic uncertainties summed in quadrature. 
A ^ jNDE of 16.28/35 = 0.47 is found using the full KEr/uq range available in the data. 
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FIG. 5: (color online) The distribution as a function of Uq^ calculated using the quadrature sum of the statistical and 
systematic errors for the hadron and photon uncertainties (a) and calculated using only the statistical errors (b). The 
calculation with an upper limit of KEr/riq < 1.0 GeV in the 20-40% centrality bin is shown with blue open circles; this is 
Range 1. The calculation with upper limits of 1.7 and 1.0 GeV in the 0-20% and 20-40% centrality bins respectively is shown 
with red * marks; this is Range 2. Horizontal lines are drawn at the location of the Xmin + 1 (solid), Xmin + 4 (dashed) and 
Xmin + 9 (dotted) for each calculation. 
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and systematic errors are used to determine the 
Figure[5](b) shows the when only statistical errors are 
included. Open circles identify the values when an up¬ 
per limit of KEr/nq < 1 GeV is applied in the 20-40% 
centrality bin. This is Range 1. It removes the re gion 
where the proton and pions deviate from riq-scaling [ 20 j . 
Another comparison, shown with * marks and referred 
to as Range 2, restricts the KEr/uq range in both cen¬ 
trality bins with upper limits of 1.7 GeV and 1.0 GeV 
for the 0-20% and 20-40% centralities respectively. This 
extends the KET/uq cut to central collisions where the 
riq-scaling is expected to remain broken When the 
KET/riq range is restricted the width of the distribu¬ 
tion increases reflecting the reduced resolving power of 
the x^ comparison when fewer data points are included. 

The optimal values for rig-scaling are located at the 
X^ minima, a value of 1.79 for all four data compar¬ 
isons. The error on the rig.^, parameter is related to the 
width of the y^ curve. It is determined from the range 
of Uq^ values where the y^ is below y^y„ -I- 1 for the Ict 
limit, y^i„ -|- 4 for the 2cr limit, and x^in + 9 the 3cr 
limit. Horizontal lines are drawn at the y^y„ -I- n values 
in Figure [5] with solid lines for the la limits, dashed lines 
for the 2cr limits and dotted lines for the 3cr limits. When 
the systematic errors are assumed to be fully correlated, 
the a = astat case, the rig’s systematic error from the 
correlation must also be obtained. The systematic error 
on the Uqj in the a = astat case is found by shifting all 
of the photon and identihed hadron V 2 values to the ex¬ 
treme maximum or minimum values in their systematic 
error ranges, re-calculating the y^ in the rig-y-space, and 
determining the Uq^ where y^ reaches a minimum value. 
The optimal nqj values and errors from this comparison 
of data points are shown with their respective x'^/NDE 

in Table [ni 


B. y^ analysis using fit to rig-scaled hadron data 

Here, a fit to the rig-scaled identified hadron data is 
used to describe the universal scaling distribution. The 
0-20% and 20-40% direct photon data is then compared 
to this function and fit using TMinuit to find the optimal 
rig-y by minimize the y^. 



FIG. 6: (color online) The 0-20% and 20-40% Au-|-Au 112 /rig 
vs KEr/riq for pions, kaons and protons are fit with a prob¬ 
ability density distribution of a Gamma function. High pr 
protons that deviate from rig-scaling in the 20-40% centrality 
bin are excluded from the fit and are not shown [^. 


assumes that the systematic errors are fully correlated 
and the photon uncertainties are limited to the statisti¬ 
cal errors, a = astat- 

To obtain 1 ' 2 /it, the rig-scaled identified hadron data 
is fit using a scaled probability density function of the 
gamma distribution, 


G{x) = A 


((x 


WW) 


( 3 ) 


where x is KEp/riq, 7 is the shape parameter, /i is the 
location parameter, [3 is the scale parameter, A is an 
overall normalization scale and F ( 7 ) is the gamma dis¬ 
tribution, F (x) = Figure [ 6 ] shows the fit 

results when the 0-20% and 20-40% Au-|-Au identified 
hadron 112 /rig data are fit to Equation [S] In the 20-40% 
centrality bin, high KEp/nq protons that deviate from 
the rig-scaled pions are excluded from the fit and are not 
shown. Table |T] lists the parameters obtained from the 
fits for both centrality bins. 
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{v2^/n 
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where ^^e direct photon V 2 and fit is the fit to the 
rig-scaled identified hadron V 2 - The y^ is summed over 
the 0-20% and 20-40% centralities comparing the V 2 fit to 
the direct photon t' 2 /ng-y where Uqj is the only param¬ 
eter. Again, the y^ minimization is performed in two 
cases to address how the direct photon uncertainty, tr.^, 
relates to the direct photon systematic errors. One case 
uses the quadrature sum of the statistical and systematic 
errors for direct photons, a = astat ® cTsys- This assumes 
the systematic errors are uncorrelated. The second case 


TABLE I: Table of the results of a Gamma distribution fit to 
the Au+Au U 2 /iig vs KEp/riq. 


Parameters 

0-20% 

20-40% 

7 

1.86 

1.62 


0.08 

0.11 

P 

1.34 

1.76 

A 

0.166 

0.34 


A TMinuit fit is used to determine the ?ig.y where the 
y^ from Equation[2]reaches its minimum value. This fit is 
performed over two ranges. Range I removes the region 
where the proton breaks the rig-scaling by applying 
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an upper limit at KET/riq < 1 GeV in the 20-40% cen¬ 
trality bin. Range 2 is restricts the KET/uq range in 
both centrality bins with upper limits of 1.7 GeV and 
1.0 GeV for the 0-20% and 20-40% centralities respec¬ 
tively. This removes the region in the 0-20% bin where 
rig-scaling is expected to be broken [2^. TMinuit finds 
the optimal Uq-y value with statistical errors. When the 
direct photon systematic errors are assumed to be fully 
correlated, the a = Gstat case, the Uq-y's systematic er¬ 
rors from this correlation must also be determined. This 
is done by shifting the direct photon V 2 values to the 
extreme maximum and minimum of the systematic er¬ 
ror range and re-fitting with TMinuit to find nq^ at the 
minimum value. The resulting Uqy values and errors 
from the TMinuit fits are shown in Table m with their 
respective /NDE. 

The low x^/NDF values under the = (Jstat 0 cTsys 
heading reflect the over-estimation of the photon and 
hadron uncertainties when uncorrelated systematic er¬ 
rors are assumed. Under the = CFstat heading, when 
only the statistical errors are used in the x^ determina¬ 
tion, the corresponding x^/^DE values are above one, 
a consequence of the underestimation of the uncertainty 
when the systematic errors are assumed to be fully corre¬ 
lated. The separation of the systematic errors into errors 
that are point-to-point independent and those that are 
correlated is needed to fully interpret the x'^/NDF val¬ 
ues in these comparisons. 

The hypothesized value of Uqy = 2 is within the sys¬ 
tematic uncertainty region when the Uq^ is determined 
from the data with (jy = Ustat in both Range 1 and Range 
2. The Uqy = 2 condition is inside of the Icr limit for the 
cr.y = Gstat 0 cfgys, Range 2 data comparison and within 
the 2(7 limit for the cr.y = fJstat 0 o'sys. Range 1 data com¬ 
parison. The Uqy values from the comparison to the fit of 
the rig-scaled hadron data are very similar to the direct 
data comparison results. A Uqy value close to 1.8 is found 
over both ranges when cr = agtat is assumed and in Range 
2 when a = Ustat^cTsys is assumed. Only the TMinuit fit 
over Range 1 produces a Uqy value that differs from 1.8, 
however, it is within 2a of the Uqy = 2 hypothesis. Of the 
eight nqy searches presented here, six are consistent with 
rig-y = 2 within Icr. The remaining two Uqy searches are 
consistent with the Uqj = 2 hypothesis at the 2a level. 
These two comparisons both use the larger KET/uq re¬ 
gion in the 0-20% centrality and a = astat 0 cfsys ■ These 
comparisons are affected by the difference between the 
pion V 2 and direct photon V 2 at KEx/nq > 1.7 GeV in 
the 0-20% centrality bin, seen in Figure[31 This difference 
between the pion and direct photon V 2 at high EEx/nq 
is the result of the increased direct photon contributions 
from hard scattering at high pr, pt > 5-5 GeV [^. 

The large systematic errors in the direct photon data 
dominate the uncertainty in the Uqy determination. Re¬ 
duced systematic errors on the direct photon V 2 mea¬ 
surement and separating the systematic errors into er¬ 
rors that are point-to-point independent and those that 
are correlated would reduce the uncertainty and improve 


the calculation of the x^ in these comparisons. Proton 
V 2 measurements that extend out to higher pr in the 0- 
20 % centrality bin, and direct photon V 2 measurements 
in additional centrality bins and collision systems would 
provide additional points for comparison benefiting this 
analysis by reducing the width of the distribution and 
improving the resolving power of the Uqy parameter. Fur¬ 
thermore, direct photon azimuthal anisotropy measure¬ 
ments at higher orders, will provide an additional 
tests to this model. The model predicts that higher or¬ 
der direct photon Vn will follow the higher-order modified 
rig scaling relation, with a universal curve in Vnlriq^"^ as 
a function of KET/riq [l^, with Uqy = 2 for direct pho¬ 
tons. 

Seven out of the eight x^ comparisons shown here find 
an optimum Uqy value of approximately 1.8. In six cases, 
the nqy = 2 condition is within \a of the optimum value. 
In the remaining two cases, the nqy = 2 condition is 
within 2a of the optimum value. These two cases are 
biased by the hard scattering contributions at high pT- 
These results, in conjunction with the similarity in the 
data seen in Figure |31 indicate that the direct photon 
V 2 data are consistent with the hypothesis of Uqy = 2 
required by the q-q annihilation production mechanism. 


III. SIMULATING THE DIRECT PHOTON V2 

To further develop the ansatz of photon production at 
confinement from coalescence-like quark-anti-quark an¬ 
nihilation, a data-driven Monte Garlo simulation is de¬ 
veloped. The crux of the direct photon puzzle is to rec¬ 
oncile the pt spectral shape with the large azimuthal 
anisotropy. In Section IIII A1 the q-q photon px spectral 
shape and V 2 are simulated with a Monte Carlo simu¬ 
lation. Rather than calculating the yields, a fit to the 
measured px distribution is performed in Section IIII Bl 
to determine if the q-q photon px shape from the Monte 
Carlo is able to describe the large excess above the Taa- 
scaled p-\-p yield seen in the data. Then the direct pho¬ 
ton V 2 is calculated by weighting the q-q photon V 2 by 
the relative contribution of the q-q photon component to 
the total direct photon yield; the T^A-scaled p-\-p contri¬ 
bution is assumed to be azimuthally isotropic. 


A. Monte Carlo of coalescence-like q-q photon V 2 
production 

The Monte Carlo consists of randomly sampling quark 
Tot values from a thermal Blast Wave distribution. The 
quark flow is implemented by calculating the quark V 2 
from a fit of the measured iig-scaled identified hadron V 2 
and then sampling the quark (/) from the U 2 -modulated 
(j) distribution. This process is repeated for three quarks 
and then co-moving requirements are applied. 

The quark’s mx is randomly sampled from a thermal 
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TABLE II: Table of optimal Uq-y values and errors with /NDF 



CT 'y 

— ^ stat 0 ^sys 

<^7 — stat 



Uq-y ± (stat) 

X^/NDF 

Uq-y ± (stat) ± {sys) 

X^INDF 

Data, Range 1 

1 7 Q + O.O 8 
‘ ^-0.27 

4.85/20 = 0.24 

1 '^Q-j-0.002+0.67 
-L- <y_0.01-0.72 

101.6/20 = 5.1 

Data, Range 2 

1.79 ±0.27 

4.53/17 = 0.27 

-I 1^0+0.002+1.09 

1- /y_o.oi_0.72 

99.5/17 = 5.9 

Fit, Range 1 

1.59 ±0.22 

3.51/13 = 0.26 

1.79 ± 0.02t(; j| 

44.67/14 = 3.19 

Fit, Range 2 

1.83 ± 0.44 

1.55/5 = 0.31 

1.88 ± 0.07to 7 ? 

34.14/6 = 5.68 


Blast Wave distribution, 

amTclyaq) 

( pt sinh(p) cos((j)) — ttit cosh(p) coshly) 
X exp I -—- 


where T is the temperature, mr = y Pt the 

transverse mass, p = tanh“^(/3s(r/i?)“) is the boost an¬ 
gle, and (j) is the azimuthal angle with respect to the 
reaction plane [2^. Further, /3s is the surface velocity, 
R is the maximum radius in the region and niq is the 
quark mass. A /3s value of 0.75 is assumed and is consis¬ 
tent with (/3) = 0.5 with a set to one. A quark mass of 
300 MeV, temperature of 106 MeV and maximum radius 
of 8.5 fm is used. The parameters of the Blast Wave dis¬ 
tribution are taken from Refs, and [2^. These Blast 
Wave parameters characterize the ttit distribution of the 
late-stage medium and therefore identical parameters are 
used for the Au-|-Au 0 — 20% and 20 — 40% centrality bins. 
The r^, y and (p values that determine the Blast Wave 
distribution are each chosen from flat distributions; r and 
y are the quark’s radius and rapidity respectively. The 
quark’s y is chosen from ±0.50 and a ±0.35 rapidity cut 
is applied to the resulting photons. The random choice 
of (j) ensures that each of the successive Blast Wave dis¬ 
tributions sample the full variation in azimuth. 

Rather than using this (p for the quark’s p, the ther¬ 
mal quark’s p is chosen from an data-driven procedure to 
reduce the simulation’s dependence on free parameters. 
This is done by using the rriT obtained from the Blast 
Wave to calculate the quark azimuthal anisotropy from 
a fit to the measured n^-scaled t >2 of identified hadrons 
shown in Figure HI Once the quark’s V 2 , V 2 q, is calcu¬ 
lated it is used to generate a 1 ± 2 v 2 q cos{2p) probability 
distribution to randomly select the quark’s p. The V 2 q is 
calculated using a fit to the measured n^-scaled identified 
hadron V 2 - A scaled probability density function of the 
gamma distribution. Equation |3l is fit to the n^-scaled 
identified hadron V 2 data as described in Section ing 

This method effectively averages the p variation within 
the Blast Wave distribution while still including radial 
boost effects. By choosing the p from the 1 ± 2v2Cos(2p) 
distribution, the measured identified hadron V 2 lnq is 
used to guide the modeled quark’s azimuthal anisotropy. 
This empirical approach to describe the quark’s az¬ 


imuthal anisotropy keeps the number of free parameters 
in the model to a minimum. One downside of this ap¬ 
proach is that the V 2 q from the fit relies on the pion data 
at high REx/nq which has increasing contributions from 
non-thermal quarks either from hard processes and frag¬ 
mentation or hard thermal coalescence [2^. This may 
underestimate the amount of quark flow at high KEt jnq. 

The random determination of the quark’s tot and p 
is repeated for the second and third quarks within the 
Monte Carlo event. The same rapidity and radius is 
assumed for subsequent quarks, and therefore the same 
Blast Wave distribution. However, a new tot value is 
sampled, V 2 q is calculated and p is sampled using the 
1 ± 2v2qCos{2p) distribution. The following co-moving 
requirements, motivated by [T^, are applied to all three 
quarks to produce a baryon and to the first and second 
quarks to produce a meson. 

Mesons: \pi — P 2 I < 2Ap, |a;i — a; 2 | < Aa; 

Baryons: |pi -P 2 I < a/2Ap, |xi - X 2 I < a/2Ax, 

|pi ±T2 — 2p3| < a/GAp, |xi ±X2 — 2x3] < v^Ax 

where Pi and x^ are the three dimensional momentum 
and position vectors of the various quarks, and Ap and 
Ax are 0.2 GeV/c and 0.85 fm respectively [l^. Quarks 
and anti-quarks that annihilate to produce photons must 
satisfy the same co-moving requirements as mesons. The 
four-momenta of quark pairs and triplets that satisfy 
the co-moving requirements are summed to create pions, 
photons and protons respectively. The hadrons and pho¬ 
tons are brought on mass shell while maintaining kinetic 
energy conservation. Figure [8] shows the amount of en¬ 
ergy taken up by the gluon to bring the photon on mass 
shell as a function of the direct photon’s KEt for the 
0-20% (left) and 20-40% (right) simulations. The z-axis 
is the number of counts and is shown with a logarithmic 
color scale. The gluon’s energy contribution is defined 
as Ej — Eqi — Eq 2 and has a value of approximately 
—600 MeV. At photon KEt < 2 GeV, Egiuon extends 
to lower energies of —770 MeV, however, the majority of 
the contribution is located at —600 MeV for all photon 
KEt values. This negative value means that the gluon 
removes some of the energy from the quarks and passes 
it to the medium when the photon is produced. Addi¬ 
tional simulations maintaining momentum conservation 
and energy conservation are also performed, however, ki¬ 
netic energy conservation best reproduces the n^-scaling 














FIG. 7: (color online) The V 2 for pions, photons, protons and 
thrown quarks simulated using the fast Monte Carlo method. 
Plots (a) and (b) are the V 2 vs KEt for the 0-20% and 20- 
40% respectively. Plots (c) and (d) are the rig-scaled results 
for 0-20% and 20-40%. The 0-20% V 2 values are scaled by 1.6 
to make the y-axis scales consistent. 


seen in the pion and proton V 2 data. Figure [7] shows 
the V 2 for the thrown quarks and simulated pions, pro¬ 
tons and photons in 0-20% (a) and 20-40% (b) centrality 
bins. The V 2 /nq vs KET/riq, Figures!?] (c) and (d), show 
that the rig-scaling is well reproduced in the simulation. 
Table IIIII displays the inverse slopes of the Monte Carlo 
Pt spectral shape when fit to an exponential in different 
Pt ranges. These are consistent with the inverse slopes 
obtained from fits to the Au-|-Au data over similar pT 
ranges [1, |1|. 


TABLE III: Table of the inverse slope of the direct photon pr 
spectral shape in different centralities and pr ranges. 


Centrality 

Pt range 

Monte Carlo 

An-f Au data [2, ^ 

0 -20% 

0.6-2.0 GeV/c 

233 ±6 

239 ± 29 ± 7 

0 -20% 

1.0-2.2 GeV/c 

251 ± 8 

221 ±19 ±19 

20-40% 

0.6-2.0 GeV/c 

233 ±8 

260 ± 33 ± 8 

20-40% 

1.0-2.2 GeV/c 

251 ± 10 

217 ± 18 ± 16 


B. Determining the yield of the q-q photon 
component 

To find the total direct photon production, a two- 
component model consisting of the q-q photon contribu¬ 


tion and the T^^-scaled p-\-p contribution is used. While 
additional photon sources are expected, these are as¬ 
sumed to be negligible compared to the q-q and Taa- 
scaled p-\-p components. The simulated q-q photon con¬ 
tributions are normalized to the measured direct pho¬ 
ton yields. The normalization constant of the q-q pho¬ 
ton component is determined from a fit to the measured 
Au-f An d, d, [13 and T^^-scaled p-\-p data d [H, |13 
using TMinuit. The normalization constant is the only 
parameter of the fit. The is calculated using the sta¬ 
tistical errors from the Monte Carlo simulation and the 
statistical and systematic errors from the data summed 
in quadrature. At low px where p-\-p reference data is 
scarce, the p-\-p yield is extrapolated from the power law 
fit obtained from d|. The normalization error on the 
q-q photon component and the systematic error of the 
T^^-scaled p-\-p ht result in systematic error band on 
the simulation. 

Figure [9] shows the resulting pp distributions for 0- 
20% and 20-40% Au-|-Au collisions. The various Au+Au 
measurements are shown in red circular symbols and the 
T^^-scaled p-\-p measurements are shown in blue square 
and cross symbols. The p-\-p ht is shown with a grey 
band, the normalized q-q photon contribution is shown 
with a purple band and the total simulated yield is shown 
with a cyan band. The error on the yield determination 
results in a systematic band on the q-q photon contri¬ 
bution which is propagated to the total simulated yield. 
Below the main hgures the ratio of the Au-f An data to 
the simulation result is shown. This ratio is ht to a hat 
line and found to be consistent with one for both central¬ 
ities, a value of 0.951 ±0.051 for 0-20% and 1.038 ±0.065 
for 20-40%. The /NDF values for these hat line hts 
are 22.8/26 = 0.877 and 32.5/26 = 1.25 for the 0-20% 
and 20-40% ratios respectively. This conhrms that the 
photons generated by the gluon-mediated annihilation of 
radially boosted quarks are able to describe the shape 
of the direct photon pp spectra for both the 0-20% and 
20-40% centrality bins. 

The total direct photon V 2 is the weighted average of 
each component’s V 2 - The r^A-scaled p-\-p contribution 
is assumed to have no reaction plane dependence and, 
therefore, a V 2 of zero. By weighting the simulated q-q 
photon V 2 by the relative contributions of the q-q photon 
yield to the total simulated yield, the total low pp pho¬ 
ton V 2 for each centrality can be calculated. Figure [TUI 
compares the simulated direct photon V 2 to the measured 
Au±Au V 2 (solid blue circles) [1|. The open red circles 
are the unweighted q-q photon V 2 generated in the Monte 
Carlo. The small black squares are the total direct pho¬ 
ton V 2 assuming uniform azimuthal production from the 
TAA-scaled p-\-p source. The relative contribution of the 
q-q photon component to the yield is shown below the V 2 
plots; this is the weight used to calculated the total sim¬ 
ulated V 2 - The error in the q-q yield normalization lead 
to the systematic error in this q-q Monte Carlo weight. 
The systematic error in the modeled V 2 is calculated from 
the quadrature sum of this normalization error and the 
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E gluon (GeV) E gluon (GeV) 

FIG. 8: (color online) The energy taken by the gluon as a function of the direct photon’s KEt for the 0-20% (a) and 20-40% (b) 
simulations. The z-axis is the number of counts and is shown with a logarithmic color scale. 
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FIG. 9: (color online) The direct photon yield versus pr for the 0-20% (a) and 20-40% (b) Au-|-Au data (red circles and 
asterisks) UMM are shown on a log-scale. The T4A-scaled p+p yields (blue squares and crosses) 0, [^ . are also shown 
including a power law fit to the p+p data (grey band) [^. The Monte Carlo yield from quark anti-quark annihilation (purple 
band) are fit to the data and are shown with the total fit yield (cyan band) found by summing the Monte Carlo yield and the 
TAA-scaled p+p fit. The ratio of the Au-|-Au data over the total fit yield is shown in the lower plots. The thick black line is a 
flat line fit to this ratio with a value of 0.951 ± 0.051 and 1.038 ± 0.065 for the 0-20% and 20-40% ratios respectively. 
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systematic error on the fit to the rig-scaled V 2 of identi¬ 
fied hadrons, with relative error values of 10% and 7% 
in 0-20% and 20-40% respectively. The model simula¬ 
tion of the total direct photon V 2 extends out to a pr 
of 3.6 GeV/c in 0-20% and 3.2 GeV in 20-40%, above 
which the simulation lacks sufficient statistics. For the 
0 -20% centrality the total direct photon V 2 agrees with 
the measured results within error bars. However, above 
a pt of 1.4 GeV/c, the simulated V 2 is systematically at 
the bottom of the error range. In the 20-40% centrality 
comparison, the total simulated V 2 agrees with the mea¬ 
sured results for pT less than 3 GeV/c, above 3 GeV/c 
it underestimates the measured V 2 - In both centralities, 
the simulated direct photon V 2 agrees with the measured 
V 2 within errors. 

IV. CONCLUSIONS 

Photon production from gluon mediated q-q annihila¬ 
tion as the system becomes color neutral is proposed as 
a large additional source of direct photons. This would 
require direct photons follow rig-scaling with an Uq-y = 2. 
The large direct photon flow measured in Au-|-Au colli¬ 
sions at RHIG is consistent with rig-scaling when rig-y = 2. 
Furthermore, in the 20-40% comparison where the high 
Pt proton V 2 lnq is seen to split from the rig-scaled pion 
result, the direct photon n2/ng follows the same trend 
as the proton. This suggests that direct photons and 
protons may experience similar transitions from the re¬ 
combination dominated intermediate pt to the higher pT 
region dominated by hard processes, comparisons of 
the direct photon and identified hadron V 2 in KET/riq re¬ 
gions where rig-scaling is seen in identified hadron data, 
find that the direct photon V 2 optimally agrees with the 
uniform rig-scaled curve when Uq^ is near a value of 1.8. 
Six out of eight of comparisons are consistent with 
Uqj = 2, the remaining two comparisons are consistent 
at the 2cr level. The two remaining comparisons include 
the 0-20% high KEr/ng region where deviations from 
rig-scaling are expected. The comparisons would ben¬ 
efit from reduced systematic errors on the direct photon 
V 2 measurement and the separation of the systematic er¬ 
rors into uncorrelated and correlated errors. Direct pho¬ 
ton and identified hadron V 2 measurements in additional 
centralities, ^Jsnn and collisions systems as well as pro¬ 
ton V 2 measurements that extend out to higher pT would 
provide further points of comparison and thus improve 
this analysis. 

A Monte Carlo simulation generates the q-q annihi¬ 
lation photon component pr shape and (j) modulation 
assuming a coalescence-like framework with quarks that 


follow a Blast Wave uit distribution and a data-driven 
V 2 parametrization. The Monte Carlo is able to repro¬ 
duce the rig-scaling of pions and protons and determine 
the q-q photon V 2 and the shape of its pr distribution. 
The resulting q-q photon pT shape with the T^^-scaled 
p-\-p photon yield is able to describe the large direct pho¬ 
ton excess seen in 0-20% and 20-40% Au-|-Au collisions. 
The simulated direct photon V 2 is consistent with the 
measured V 2 in the 0-20% centrality bin but systemati¬ 
cally low. In the 20-40% comparison the simulated direct 
photon V 2 is able to reproduce the measured direct pho¬ 
ton V 2 at Pt less than 3 GeV/c but underestimates the 
V 2 at higher pT as the T^i^-scaled p-\-p contribution be¬ 
comes significant. The addition of thermal hard quark 
pairs would likely contribute to additional yield and flow 
for Pt values above 3 GeV/c [1^. Future work would 
benefit from a more robust hydrodynamic calculation of 
the flowing quarks near the phase transition with yield 
estimates. This is particularly important as the deter¬ 
mination of the quarks V 2 from the iig-scaled identified 
hadron V 2 is expected to falter at high pT as non-thermal 
production mechanisms such as thermal hard coalescence 
and fragmentation from hard interactions contribute to 
the pion yield. 

This paper has focused on the published 
= 200 GeV Au-f Au 0-20% and 20-40% di¬ 
rect photon Pt and V 2 distributions. Future work to 
simulate the q-q photon contributions in more peripheral 
collisions is promising. Additionally, the higher orders 
of the direct photon flow presents a new quantity to dis¬ 
tinguish between the different photon processes. Given 
the soft-gluon mediated q-q annihilation production 
mechanism ansatz, the Vn for direct photons is expected 
to be similar to the pion Vn at pt less than 3 GeV/c for 
higher orders of n. This model predicts that higher-order 
Vn rig-scaling laws seen with identified hadrons [T^ will 
remain valid for the direct photon Vn where the nq^ = 2. 
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